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Abstract: We study squeezed limit fwi generation by excited initial inflationary 
states in a model independent way. We restrict "excited" to mean a Bogoliubov trans- 
formation of the Bunch Davies state. We simultaneously impose the constraints that 
the observable power spectrum is nearly scale invariant over at least three decades and 
that the observable modes today be subhorizon at the beginning of inflation while not 
causing significant backreaction. We show that most excited initial inflationary states 
for single field inflationary models with negligible superhorizon evolution do not pro- 
duce an observable squeezed limit Jnl- The case in which one mode is in the Bunch 
Davies state while the other two modes are in an excited state with 0.01 < \(3k\ < 0.1 
may generate a squeezed limit /jvx which is detectable with future experiments. 
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1 Introduction 

The standard single field inflationary scenario whereby a single scalar field undergoes 
slow roll inflation has been remarkably consistent with observational data[l]. The usual 
assumption is that the fluctuations that we observe in our universe today were initially 
in a Bunch Davies state. However, recently there have been attempts to study models 
of inflation which have an initial state which is related to the standard Bunch Davies 
vacuum through a Bogoliubov transformation[2-ll]. We are restricting our usage of 
the word "excited" to mean only states which have a well-defined non-zero Bogoliubov 
coefficient 0^. We also are not considering models with superhorizon evolution such as 
those studied in [12, 13]. It was shown in [14] that pre-inflation background anisotropy 
models, in the so called non-planar limit, result in an excited initial state with rescaled 
Bogoliubov coefficients given in equation (3.9) of that paper. Therefore these models 
should be constrained by our analysis as well, the only possible difference being that 
the Bogoliubov coefficients depend on a special direction. We will show most models of 
inflation with excited initial states and negligible superhorizon evolution cannot con- 
sistently result in an observably large squeezed limit /nl- 
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In section 2 we prove that in general, based upon the observation of three decades 
of nearly scale invariant power spectrum, most excited initial states cannot be obser- 
vationally distinguished from the standard Bunch Davies state by the measurement of 
squeezed limit f^L- We will consider both the mixed case whereby today's observable 
modes are composed of both Bunch Davies modes and excited modes in addition to 
the unmixed case whereby today's observable modes are composed solely of the latter. 
In order to make our argument less abstract, we explicitly demonstrate this for the in- 
teresting case of large excitations in section 3. We will begin section 3 by reviewing an 
excited state scenario in the literature that is physically motivated by the introduction 
of an initially kinetic dominated era prior to the inflationary era [7] . We then explicitly 
show that in this model the squeezed limit Jnl for the case in which all modes under 
consideration are in the excited initial state is not observable. In section 4 we conclude. 



2 Excited Initial State f^L Undetectability 

In order to study the squeezed limit f^L for single field inflationary scenarios in a model 
independent fashion, we will simultaneously impose four constraints: 1 

(1) Modes that we observe today were subhorizon at the beginning of inflation. 

(2) Slow roll inflation occurred undisturbed by fluctuation backreaction. 

(3) Any superhorizon evolution was negligible. 

(4) There are at least three decades of nearly scale invariant observable power spectrum. 

In sections 2.1 and 2.2 we impose these constraints to show that the only scenario 
which may result in an order unity Jnl is the case in which one mode is in the Bunch 
Davies state while the two modes are in an excited state with 0.01 < \(5k\ < 0.1. In 
section 2.3 we will parameterize these constraints explicitly in order to demonstrate 
their robustness. 

2.1 Mixed Case: Bunch Davies Modes and Excited Modes 

For general excited states, the power spectrum of comoving curvature perturbations 
(71) is given by 

with | ctfc | 2 — \(3k\ 2 = 1 (see for example, [6]). The near scale invariance of the observed 
power spectrum therefore prohibits the case of mixing Bunch Davies modes with excited 

1 Also see Raphael Flauger's presentation at the workshop Critical Tests of Inflation Using Non- 
Gaussianity at http:/ /www. mpa-garching.mpg.de/~komatsu/meetings/ng20I2/. 
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modes unless \/3k\ < 0.1 [2]. The leading contribution to f^L in the mixing case for small 
is of the form e\fik\(k\/kz) [6] with (ki/k 3 ) < 10 3 , where hi is the short wavelength 
mode and ks is the long wavelength mode. This assumes that the two short wavelength 
modes are in an excited state and the third long wavelength mode is in the Bunch 
Davies state, otherwise one does not obtain an enhancement in the squeezed limit. 
Constraining e in terms of the tensor to scalar ratio, r < 0.13[1], we obtain 2 
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W + I3 k \ 2 . (2.2) 



The constraint becomes e < 10~ 2 for the mixing case of \(5k\ < 0.1, resulting in a 
maximum squeezed limit Jnl contribution, up to numerical coefficients, of order unity. 
For \/3k\ < 0.01, however, we obtain /jvx ~ lO^ 1 which is undetectably small. Therefore 
the mixed case does not result in an observationally distinguishable squeezed limit f NL 
from the standard Bunch Davies result if \p k \ < 0.01, but may for 0.01 < |/3fc| < 0.1. 

2.2 Unmixed Case: Only Excited Modes 

Now we restrict our attention to the case of all observable modes being initially in an 
excited state. That the modes are observable today dictates that they were subhorizon 
at the beginning of inflation, H < p, where p is the physical momentum. Addition- 
ally, we require that the excited state modes do not generate enough backreaction to 
spoil slow roll inflation nor bring into question the validity of perturbation theory by 
imposing that the fluctuation energy density be much less than the total energy den- 
sity, (pn) <C H 2 Mp. Computing the energy density of the fluctuations explicitly using 
adiabatic subtraction and assuming that the Bogoliubov coefficients are approximately 
scale invariant 3,4 {\Pk\ ~ |/?|) gives a leading term which is quartic in momentum, 



|2 



(pk) = e y o ^ \R k \ 2 - i/ BD , k i 2 + r) (\R k \ 2 - iw) 



(2.3) 



8tt 



2 Hcrc we have assumed that the tensor spectrum is left unchanged by setting the initial state of 
scalar perturbations to be something other than Bunch Davies. Using an initially excited state for the 
tensor perturbations would further suppress e and therefore further suppress f^L as we shall see. 

3 Except for in the case of <C 1 which is too small to have its scale dependence constrained 
by power spectrum observations, one may take any value for |/3| that well represents the decades of 
interest and it will not alter our conclusion. 

4 This differs very slightly from the analysis conducted in [4, 6] whereby it was assumed the coefficient 
would exponentially decay in momentum. However this difference will not invalidate our discussion 
of backreaction since the modes whose contribution survives adiabatic subtraction have an essentially 
constant exponential term which is near unity. 
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Here, puv is the highest physical momentum mode which is in an initially excited 
state and Rk is the comoving curvature perturbation mode function related to the 
Bunch Davies mode function /bd,1c through the standard Bogoliubov transformation. 



The upper bound on the physical momentum is then given by p <C \/10HMp/\/3\. 
Therefore combining our subhorizon lower bound with our backreaction upper bound 
we restrict the range of allowed momenta modes to be, 

H < p < ^/lOHMp/lPl. (2.4) 

The range of observable modes today for which scale invariance holds spans at 
least three decades [1, 15] and therefore we require the conservative case that there be 
at least three decades of allowed physical momenta range, 



, JlOHMp/\f3\ . 

10 3 < TT 1 — . (2.5) 

H 

However the observable amplitude of the power spectrum [1] constrains the Hubble 
parameter H in terms of the slow roll parameter e and Bogoliubov coefficients and /3k, 
(Mp/H) 2 < 10 7 e _1 |o; + /3| 2 . 5 Using this correspondence the constraint (2.5) becomes, 

10 3 |/3| 2 e< |a + /3| 2 . (2.6) 

There are two cases to consider, that of large \/3\ and that of small \/3\. For the 
case of large \/3\ > 1 the inequality restricts e to be very small, e <C 1CT 3 . The lead- 
ing contribution to squeezed limit f^L in this case is proportional to e(ki/ks)[6], with 
(&1/&3) < 10 3 . This results in f^L, up to numerical coefficients, being forced to be 
much smaller than order unity, f^L <C O (1). 

For the case of small \(3\ < 1 the constraint simplifies to |/3| 2 e <C 1CT 3 , which upon 
substitution of (2.2) results in |/3| 2 <^ lO^ 1 . As in the mixed case, the largest value of 
e is chosen since we would like to maximize f^L which has a leading contribution of 
the form |/3|e(/ci//c 3 ) [6]. 6 This is clearly undetectable as well since we again find, up 
to numerical coefficients, f^L 

Therefore inflationary models with all observable modes as excited states do not 
have a squeezed limit f^L which is observationally distinguishable from the Bunch 
Davies initial state. We will now parameterize the inequalities used in the proof in 
order to study its robustness. 



= 2^ Pn{k) = 8^M^ |afe +/?fe|2 = 2A X 1(r9 ' Therefore [-Jf] ~ 5 x 10 6 |a fc +/3 fc | 2 A- 
6 Up to numerical coefficients, |/3| 2 e < fO -3 may be written as f NL < |/3 1 elO 3 < y/elO 3 ^ 2 , which 
motivates the choice of choosing the largest e. 
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2.3 Unmixed Case: General Parameterization 

In order to develop an intuition for the robustness of the proof previously given, we 
now repeat the argument directly parameterizing the inequalities. In general one may 
take, 

(1) Subhorizon constraint: A SU ^H < p. 

(2) Backreaction constraint: (pn) < A^^^Mp. 

(3) ndec nearly scale invariant observable decades. 

This parameterization leads to the constraint which is analogous to (2.6) of 

10 4 " d -" 9 |/3| 2 e < ^\a + (3\ 2 . (2.7) 

^sub 

For the case of large \fi k \ > 1 the constraint (2.7) simplifies to e < 10 9 " 4ndcc (A hack /A^ nh ) . 
Therefore the leading contribution to f^L is given, up to numerical coefficients, by 

Inl ~ e(h/h) < 10 12 - 4m - (A back /A 4 ub ) . > 1) (2.8) 

For the case of small \p k \ < 1 the constraint (2.7) simplifies to \(3\ 2 < lO 11-4 ^- (A hack /A 4 suh ) , 
where we have once again used (2.2) to constrain e < 1CT 2 and chosen the largest value 
of e in order to maximize f^L- Up to numerical coefficients, /atl is then given by 

f NL ~ \P\e(ki/h) < ^10i3-4n, cc (A hack /At h ). (\(3\ < 1) (2.9) 

In order to gain some intuition for these numbers, consider the values given in 
Table 1. The chosen values for the parameterization variables are very conservative, 
physically corresponding to modes being barely inside the horizon at the beginning of 
inflation and the energy density of fluctuations being as large as one tenth the total 
energy density. For the number of nearly scale invariant observable decades, we note 
that three decades is a conservative lower bound [15] and therefore demonstrate the 
dependence as we allow more scales to be considered. The values found for f NL are 
always less than order unity in both the case of large \(3\ (2.8) and small \(3\ (2.9). 



-^■sub ^back ^dec 


Large \/3\ f NL (2.8) 


Small \(3\ f NL (2.9) 


2 0.1 3 
2 0.1 3.5 
2 0.1 4 


< o(io- 3 ) 

< c(io- 5 ) 

< o(io- 7 ) 


<0{\Q- 1 ) 

< 0(1O- 2 ) 

< o(io- 3 ) 



TABLE 1: Typical f^L values for conservative parameterization variables. 

Clearly, for any realistic model of inflation one will not be able to generate an observable 
squeezed limit f^L for models in which all observable modes today are in the excited 
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state. It should be emphasized that our conclusion of an undetectable squeezed limit 
/ml for the unmixed case required only the four constraints discussed at the beginning 
of section 2. Therefore the analysis is very model independent. However, as a concrete 
example we will present in the next section a specific model with physically motivated 
large Bogoliubov coefficients. 



3 Example: Case of Slow Roll Enhanced \(3k\ 

3.1 Review of the Scenario 

We now explicitly show an example to illustrate that inflationary models with excited 
state vacua that have slow roll enhanced Bogoliubov coefficients cannot have observable 
consequences today. These large Bogoliubov coefficients are motivated by allowing for 
a kinetic dominated era before inflation begins. Following the model in [7], expansion 
of the universe is driven by a single scalar field (the inflaton) minimally coupled to 
gravity. The inflationary potential is of the form, 

f , < fa (Kinetic Era) 

W \ A 4 (l - a(p) , (pi < (p < (p f (Inflationary Era) ' 1 ' ' 

where a<p <C 1 for fa < (p < (pf. We will assume that post- inflationary physics does not 
affect our results for the entirety of this study. 



We are interested in the behavior of the comoving curvature perturbation at the 
transition point fa. In this study we are focused on modes which will potentially be 
observable today and therefore restrict our attention to classifying modes which are 
subhorizon at the onset of inflation. Depending on how the modes experience the 
transition they can be divided into four regimes according to the physical momentum 
scale, p, and characteristic transition time, r: 

H <C p <C r _1 : (Sudden Transition) 
H <C p ~ t -1 : (Interpolating Regime) 
H <^ r^ 1 <^ p : (Adiabatic Transition) 
Mp < p : (No Transition) 

The mode function for the kinetic 7 and inflationary periods for modes that undergo 
a sudden transition was reported in [7], 



K h (t) = 




, (Kinetic Era), 



7 The kinetic era mode function is consistent with the form of the Mukhanov variable used in [16] 
motivated by taking the Bunch Davies vacuum for very short wavelengths. 
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K k {rj) 



1 H 



[a(k) (1 + ikrj) e~ ikv + f3(k) (1 - ik V ) e ifc "] , ( Inflationary Era) 



/ Sudden 
\ Transition 



a 



(^) ~ V^Le 1 



(2) 



3(1 +iw)H{ 



( 2 ) ( 

2 



— > 



(3.2) 



//.f ('.';) + 3 (1 - iw) h[ 2) 1 " 



Here, H is the Hubble parameter during inflation, a is the scale factor with <2j as its 
value at the beginning of inflation, w = k/aiH, 77 = —1/aH, and e = 2 /2H 2 Mp. The 

(2) (2) 

Hankel functions are denoted by Hq and f/j . In the subhorizon limit, w ^> 1, these 
expressions reduce to, 

31 

4? 



|a(fc)| 2 = |0(*Of 



Arg 



2w. 



(3.3) 



3.2 Computing / 



NL 



For the sake of clarity, the explicit computation of J'nl in the standard squeezed limit 
where we have required all three modes satisfy (3.3) is located in the appendix. All 
observable modes in this case must be in the excited state since otherwise the large 
would violate the observed power spectrum scale invariance. Since this is a model 
with large \(3k\, the result of that computation (5.5) is an f NL that is proportional to 
the slow roll parameter e as discussed in section 2, 



NL 



24 



(3.4) 



In order to satisfy (2.6) we require that e 1CT 3 and therefore at most we have 
Inl <S 10 , which is unobservable. Therefore this model cannot provide a detectable 
squeezed limit f^L an d therefore cannot be experimentally discriminated from the 
standard Bunch Davies scenario. 8 



4 Conclusions 

We have argued that in order to be consistent with the nearly scale invariant observable 
power spectrum, one must require that all decades of observable scale today must either 

8 The tensor spectrum would not be affected by the transition to leading order since the total energy 
density is constant across the transition and therefore would not provide any additional evidence in 
favor of the pre-inflationary kinetic era. 
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be an excited initial state, or a combination of the Bunch Davies and excited initial 
states with small \(3k\- The mixed case of allowing for one observable mode to be Bunch 
Davies and two of the modes to be excited results in a small contribution to Jnl for 
\Pk\ < 0.01 while the contribution for 0.01 < \(3k\ < 0.1, however, may indeed be order 
unity. Allowing for the three decades to be in the same initial excited state imposes 
strong constraints on f^L once both the constraints of subhorizon scale and negligible 
backreaction are taken into account. We find in particular that the constraints for 
small |/3 1 models result in a heavily suppressed \(3\e and for large \/3\ models a heavily 
suppressed e. The leading contribution to the squeezed limit /nl for each respective 
case is proportional to these terms, resulting in a prediction that is undetectably small. 

Therefore most models of single field inflation with excited initial states and negli- 
gible superhorizon evolution cannot consistently produce an observably large squeezed 
limit Jnl for observable modes. It is only within the small region of parameter space 
whereby two modes of the squeezed limit are in an excited state with 0.01 < \(3k\ < 0.1 
and the third is in the Bunch Davies state that one may obtain an /jvl of order unity. 
If the constraints on the scalar spectral index or the tensor to scalar ratio improve with 
future data, the constraints on and e will respectively become more restrictive and 
may allow us to rule out detecting non-gaussianity from all excited state inflationary 
models with negligible superhorizon evolution. 

NOTE: As we were completing this paper, R. Flauger, D. Green and R. A. Porto 
informed us of their work [17], in which they reach similar conclusions. 
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5 Appendix A: Computing /jvx for the Case of Slow Roll En- 
hanced 

To compute the bispectrum up to tree level, we follow the approach in [18]. The opera- 
tor TZ(k, t) is defined the same way as in the standard case, but the mode function (3.2) 
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is used instead of the standard Bunch-Davies mode function. The following expression 
is obtained in the limit where all the ki modes satisfy (3.3), 



(ft kl ft k2 ft k3 ) = (2tt) 3 <5 (kj + k 2 + k 



H L2 1 

3) ^ n(2fe?) 



Bk 1 k 2 k 3 
Pk 1 k 2 k 3 



36 ]1 cos 2 (w a ] 



o 4> v* ^1 



I" O (72 S COS 2 fli)0 2^1 



H<j> 4^ cos 2 (iOi,) 2 H 2 | 4^ cos 2 (u) b ) Z^ cos 2 (u> 6 ) 
b | o o^c 



^ H 2 S ^a^6 



3+S c cos(6>t-0 e )+2£ e cosf 



fci+fc 2 +fc 3 ' 1 fc i fc 2fc3 

1+2 cos 6>3-cos(6>i+e2)+cos(6>3-6>i)+cos(6> 3 -6'2) 

fcl+fc2-fc3 



+ cyclic 



Here, 9 b = 2w b = 2k b /(aiH) and e t = 6i + 02 + #3 



(5-1) 



The relative magnitude of the bispectrum in comparison to the power spectrum 
can be written in terms of /tvl [19] as, 



B n (h, k 2 , k 3 ) = x [P n {h)P n (k 2 ) + P n (k 2 )Pn{h) + P n (h) P n (ki)\ 

5 



(5.2) 



For our case, the power spectrum for modes satisfying (3.3) has been derived in [7] to 
be, 



Pn{k) 



i m 2 

2k 3 2M 2 e 2 



cos 2 (w) . 



(5.3) 



In the squeezed limit we use the relation, 



12 



lim B n (ki, k 2 , k 3 ) = —fNLpR.(ki)Pn(h), 

fci~fc 2 >fc 3 5 



(5.4) 



to obtain the squeezed-limit /jvx- This gives us 
fNL(knh) 



Mt ||i+/I 



24 \k 



(5.5) 



Here, the term / is a factor proportional to the sum of cosines of various combinations 
of w a = k a /(a,iH), which depends sharply on the choice of modes (k) and can be taken 
to average out in measurements to zero. 
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